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Abstract 

The isoscaler giant monopole resonances (ISGMR) are computed using the canonical-basis time- 
dependent Hartree-Fock-Bogoliubov theory (Cb-TDHFB) with five kinds of Skyrme parameter 
sets (SGII, SkM*, SLy4, SkT3 and SkI3). To extract the nuclear matter property from finite 
system, ISGMRs of N=Z (Z =20 - 50), isobar even-even nuclide for A=100, 132 and Sn isotopes 
are analysed systematically. The magnitude relation of nuclear incompressibility-parameter (K 0Q ) 
among Skyrme parameter sets, can be corresponded to the peak positions of GMR in spherical 
isotopes over A=80. The parameters (AT sur f, K T and Acoui) which appear in expansion of the finite 
nucleus incompressibility Ka, are determined for each Skyrme parameter. From the comparison 
experimental data whole mass region and the present results, they indicate that the isospin de¬ 
pendent term K r is filtered as -305T10 MeV. The incompressibility parameters of infinite system 
corresponding to our results are iF(?°=-340±35, K oa = 225±11, and A" sym =-138±18 MeV. 
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I. INTRODUCTION 


To extract an equation of state (EOS) for nuclear matter from finite nuclear system, is 
one of most important task given to nuclear physics. The EOS is a very important topic 
to connect the nuclear physics to the astrophysical objects such as nucleosynthesis, neutron 
star, and so on. The EOS is often expressed in the expansion around symmetric matter, as 
follows. 
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where 5 = (p n — p p )/p is an asymmetric parameter which separates a symmetric matter 
(SM) EOS £ and symmetry energy S. They are expressed in g = (p — po)/3po, which is an 
expansion around the nuclear saturation density p 0 at which the pressure of nuclear matter 
is zero. 

When the EOS is expanded as Eq. (JTJ), there appear characteristic parameters: the binding 
energy per nucleon £ 0 , the incompressibility of SM K^, the skewness parameter of SM Q 0 , 
the symmetry energy J, the slope parameter L and the symmetry incompressibility K sym , 
and so on. The £q and J are respectively equal to the £ and S at po- The other parameters 
can be obtained from the density derivation of £ and S as follows. 
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Thus the EOS parameters can be easily obtained at when the energy density functional 
(EDF) E/A[p\ is chosen. Nucleus is a too much small system to extrapolate directly to the 


infinite nuclear matter 


Q. 


Therefore, we evaluate the EOS with the help of effective inter¬ 


action (Skyrme, Gogny, relativistic mean field) 


-9j. However, it is not easy to determine 
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acceptable values in nuclear structure and in astrophysics, at same time. Currently also 


still, a study of the 
also experimentally 


OS using effective interactions has been progressed theoretically and 


10 


13] . The uncertainties of £ 0 and J are expected to be small: £ ( 


16 MeV which has appeared in the Bethe-Weizacker mass formula, and J=32±3 MeV 
The does not have so large uncertainty due to the consistency among experiments and 
theoretical prediction: A' oo =230±30 MeV |6j. The slope parameter L which will be strongly 


related to nuclear dipole mode, has been well studied from 


instance the relation among the neutron-skin thickness (l3. 1 1 5H 1 8t l , py gmy dipole resonance 


19] and polarizability 
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re many points of view, for 
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20], although it has been 


[ltil ]. giant dipole resonance (GDR) 

yet floated: A=58±18 MeV ( 5 , 9|. These untiring studies narrow downs the range of the 
EOS parameters, however the K sym especially has a large range of values. Basically, it is 
difficult to connect between EOS parameters and experimental values directly, therefore we 
usually take the procedure: to search an EDF to reproduce experiments and then to extract 
the EOS parameters from the EDF. In a present work, we evaluate incompressibility pa¬ 
rameters from the isoscalar giant monopole resonance (ISGMR) of finite nuclear system. In 
order to estimate the parameters independently of the speciality of each nucleus, we use the 
finite nuclear incompressibility A^-expansion to analyse them, although we can compare 
the ISGMRs in theory and in experiments. 

If the energy of ISGMR Eqmr is represented in the root mean square radius of the nucleus 
and the GMR can be regarded as a single phonon mode, Eqmr is written in 


h 2 K, 
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where K A is defined as an incompressibility of finite nuclear system [1]. Although the K A 
can not directly equal to Koo at a limit of A —y oo, it will bring the relation between the 
parameters of finite and infinite system. The K A is expanded around K <*, as 
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When incompressibility parameters are extracted from the experiments, the expansion (J7J) 


is often used 
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14 ]. The surface and Coulomb term (AT sur f, A'coui) are estimated by 


infinite EOS parameters in previous studies, which will be confirmed how available in this 
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work. The isospin term K r has not been clear, which is a most important key to know the 
information of nuclear matter incompressibility using Eq. (17|) . 

This paper is organized as follows. In Sec HO we present the methods to calculate ISM 
mode and to evaluate the centroid energy of ISGMR E G M r. In Sec IIII1 at first, we show 
partially the strength functions of N = Z nuclide to check the behaviour of them including 
the results of deformed nuclei. After that, we extract the EOS parameters of finite nuclear 
system by the chi-square fitting with Eq.([7|), from the Ka evaluated with Eq.([6]) and E G mr. 
In particular, we investigate N = Z nuclide for K s ur f and A'coui, and isobar nuclide for K r . 
Furthermore, to confirm the expansion itself and with our coefficients, we compare them with 
the experimental E GMR and calculated E G mr for Sn isotopes. In Sec lIVl we also compare 
the experimental measurements for a whole mass region with A=24—238, in order to narrow 
down the candidates of effective interaction. And we mention the mass dependence of E G mr 
in experiments. The present work is summarized in SecfVl 


II. FORMULATION 


To access the GMR, we apply the Canonical-basis time-dependent Hartree-Fock- 


Bogoliubov (Cb-TDHFB) theory 2l| in three-dimensional (3D) coordinate space w’ 


be successfully applied to the study of the dipole 


0,G] 


rich can 


and quadrupole 


23 


24 | modes 


of many isotopes, systematically. The Cb-TDHFB can describe the dynamical effects of 
pairing correlation in fully self-consistently. The Cb-TDHFB equations are derived from the 
full TDHFB equation represented in the canonical basis {</>;(£), which diagonalize a 

density matrix, and by assuming the diagonal form of pairing functional. The Cb-TDHFB 
equations compose the time-evolution equations for the canonical pair {(fti(t), </y(f)}, its oc¬ 
cupation probability pi(t) and pair probability 


1 = ^ ~~ ( 8 ) 
i Jt Pl ^ = ~ 

i = (77!(t) + 77f(t)) Ki(t) + At(t) (2pi(t) - 1), 


where the phase of canonical basis is chosen as r/i(t) = (<f)i(t)\h(t)\<f)i(t)), and the h(t) and 
A i{t) are the single-particle Hamiltonian and the gap energy, respectively. 
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Due to the appearance of deformed ground state in our subjective nuclide, we should 
choose a flexible calculation space. We use the 3D Cartesian coordinate-space representation 
for canonical basis, a; t) = (r, a\<pi(t)) with spin a = ±1/2. The condition of calculation 
space is discretized in a square mesh of 1.0 fm inside of a sphere of radius 15 fm for all nuclide 
in present work. 

To apply our method to systematic investigation, we choose the Skyrme EDF to ph- 
channel and the simple pairing functional form to pp(hh)~ channel: A i(t) = ^2 k G k iKk{t) 
where Gki is constant in real-time evolution, as same as Ref. 2l|. Our choice of five Skyrme 


parameter sets is SG11 


25|. SkM* 26], SLy4 27], SkT3 28] and SkI3 29]. The reason to choose 


them is not only their usefulness, also corresponds to the limitation of iC oo =230±30 MeV 
indicated in Ref. js] 


A. Linear response calculation with Cb-TDHFB 


In order to induce monopole responses, we add a weak instantaneous external held 
V ext (r, t ) = £F(r)6(t) to initial states of the time evolution. Here the isoscalar monopole op¬ 
erator acting on nucleus is given as F/q = r 2 Yoo- The amplitude of the external held is so cho¬ 
sen to be a small number £ = 2 x 10~ 3 fm -2 to guarantee the linearity. The strength function 
S(E] E0) can be obtained through the Fourier transformation of E{t) = (T(t)|F/ ; g|'k(t)): 


S{E\ EO) = ^|{<I>„|FM|<t 0 )| 2 i(B„-r) 

n 

r°° 

= —Im / [E(t) - J’(O)] e< E+lV/2)t dt, 

tt? Jo 


(9) 


where |T 0 ) and |\l/ n ) are the ground and excited states, respectively. T is a smoothing 
parameter set to 1 MeV for whole nuclide in present. 


B. Evaluation of mean energy for GMR 

We need a procedure to compute the mean energy of GR, in common among subjective 
nuclide. To compute the centroid energy E of GR, we use mi/m 0 , although there are some 
evaluations or 
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The mi/rriQ is computed as follows in present work. 


_m 0 /> e e dE ES(E; EO) 

U'GMR = — —— e -;- 

mi f >£ e dE S(E ; EO) 


( 10 ) 


where < e and > e are an upper and under cut-off energy respectively. They should be decided 
with more carefully, because the Eq mr is sensitive to them. In this work, we decide them 


as: 


< e = Ec + 7.5 and > e = E q — 7.5 MeV, where Eq = 8071 4 / 3 . The empirical formula 


Ac is found in a droplet model [30]. Purposely we chose this way, because Eq. ([6]) to relate 
the Ka with nuclear response, which is based on the one phonon picture in other words the 
GMR is assumed as one mode. 


III. RESULT 

We evaluate the mean energy of GMR using Eq. flTOlh and from them the finite incom¬ 
pressibility Ka is also evaluated according to Eq.(l6|) with using calculated (r 2 ) in Table [T] 
To determine the expression parameters of Ka according to Eq. (]T|) . we proceed a following 
way step by step. First, we determine the K sur fandKc ou \ using the results of N = Z nuclide 
in which the isospin term K T does not contribute to Ka- Second, we determine the K r 
while using the A" surf and Acoui fixed in the first step. To obtain the K r for each Skyrme 
interaction, we analyse isobar nuclide for both A =100 and 132 at same time. Lastly, to 
confirm that the expansion of Ka with determined parameters reproduces the results of Sn 
isotopes, we calculate the centroid energies Agmr according to Eq.(JE]) with Ka and (r 2 ) in 
Table [TO and compare them with actually calculated Agmrs. 

A. N = Z 

ISM strength functions of even-even N = Z nuclide from 40 Ca to 100 Sn with SkM* are 
shown in FiglU Chain, dashed, doted and thick lines show the results of Z =20, 30, 40 
and 50, respectively. We can see the broad strength distribution of 40 Ca in heigh energy 
around 21 MeV. The distribution becomes localized and its centre shifts to low energy, as 
mass number increase. In these strength, split distributions can be seen in thin lines which 
are corresponding to 48 Cr and ' 2 Kr. The split is caused by the coupling monopole with 
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quadrupole excitations due to the deformation. Typically the quadrupole GR appear in 
lower energy than GMR, thus the -Egmr of a well deformed nucleus shifts to low energy. 

There are some strengths in vicinity of zero energy, which corresponds to numerical 
spurious mode due to the detail of mesh size and of time step. They are excluded from the 
estimation of E GMR in Eq. dTUD . 



FIG. 1: (Colour online) Strength functions of isoscalar monopole vibrational modes of N = Z 
even-even nuclide from 40 Ca to 100 Sn. 

Figure [2] shows the Eqmr for N = Z nuclide with five Skyrme parameters. Filled symbols 
means results of spherical nuclei or the nuclei which have small deformation (|/3| < 0.1), 
open ones means those in deformed nuclei. Over A = 80, the trend and relation among the 
results of each Skyrme parameter become clarified. The behaviour of Fgmr hi deformed 
nuclei diverges from the trend of spherical. The order of -ERmr can be almost corresponded 
to the order of K^ magnitude (refebdlerThb lllll) . 



N = Z 

FIG. 2: (Colour online) Mean energies of ISGMR for N = Z even-even nuclide, which are computed 


by Eg. (1101) using the strengths with five Skyrme parameter sets. 
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Figure [3] shows Ka obtained by Eq.()6|). Same as Figj2} the hied and open symbols 
correspond to the results of spherical and deformed nuclei, but we use the star symbol for 
double magic (DM) nuclei ( 40 Ca, 56 Ni, 100 Sn). The behaviour of Ka well corresponds to 
that of Eq mr, thus the trend and the relation among interactions are stable over A = 80, 
and the results of deformed nuclei clearly have difference from spherical nuclei. 



FIG. 3: (Colour online) Finite nuclear incompressibility Ka for N = Z even-even nuclide. They 
are estimated with Eq.©. 


To obtain the expansion coefficients (iF su rf, K T , Rc ou i), we analyse our results according 
to the Ka expansion in Eq.©. The results of deformed nucleus are excluded from our 
analysis, because they have clearly a different trend from those of spherical nuclide. If they 
can be included into the analysis, we will need the way to separate quadrupole and monopole 
modes. Here, two cases are considered: (i) excluding DM nuclei and (ii) including them. 
Our purpose is to extract nuclear matter properties from nucleus. The effects due to the 
special nuclear structure such as a modes coupling in deformed nuclei, should be excluded 
from the matter property analysis. As mentioned in Ref.[32j, the magicity effects in the 
incompressibility will appear, which should be confirmed in the comparison the (i) and (ii). 

In this section, we fix the K SUI f and Ac ou i expansion coefficients which are listed in the 
Table ED The surface term K suri is often estimated as an opposite sign of K^ J 2 , gj, 
and the Coulomb term A'coui is estimated in Ref.jl, 2] as 
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Our | AT sur f | are larger than \K 00 \, which have been mentioned already in Ref. 


( 11 ) 

Although an 


assumption AT surf = — K^ is sometimes used for non-relativistic models in previous analyses 















, it is not suitable actually and the difference over 30% from it might cause 
a serious missing in nuclear property. In both (i) and (ii), Ac ou i closes to the Acoui- AT sur f 
and Acoui in (ii) are a little weaker and stronger than those of (i), respectively in most 
interaction-cases. The effect of DM on Ka is regarded as small as shown in FigJ3} although 
small kinks appear at 56 Ni. The nuclear magicity is not so sensitive to the A% excluding 
light DM nuclei. 
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B. A = 100,132 


In this section, we determine the isospin term K r from isobar nuclide A=100 and 132, 
while using A" sur f and Ac ou i fixed in previous section. In same as Sec. IIIIA1 we exclude 
the deformed nuclei from the analysis. Figure ® shows the Ka of the spherical isobars with 
A=100 and 132 with respect to isospin asymmetry (N — Z)/A, in which the vertical chain 
line separates A=100 and 132. The root mean square radii and quadrupole deformations of 
the isobars are listed in Table [H 

The Ka in Figj4]has a parabolic shape in (N — Z)/A which corresponds to the expansion 
Eq.flU), however the centre of the parabolic function is not always at N=Z. In this analysis, 
we also consider the two cases: (i) with DM and (ii) without DM. The isospin term K T 
obtained in the cases are listed in Table ELD The effects of DM are not so large also in K r 
excluding the results of SkI3. When we exclude the result of DM from SkI3 results, the 
points which can be used in the analysis are only four, therefore the analysis ambiguity 
becomes large. To determine the A%, the number of isobar nuclide is essential. 

In several papers, an isospin dependence of incompressibility for nuclear matter is es¬ 


timated at saturation density p 0 with a small isospin asymmetry 
written in 


14|, which can be 


A;°° = K sym -6L- ^LL. (12) 

-^oo 

The K°° can not be regarded as the finite incompressibility K T , which is mentioned also 
in Refill), although the strong correlation between K r and K^° can be expected naively. 
The K^° is often used to expand the K r : K T = K™ + A'® urf A” 1 / 3 . The A'® urf could not be 
estimated because the number of isobar-chain sample is only two. The absolute value of K r 
is usual smaller than K^° in present work and the correlation among them does not seem 
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simple. 



0 0.05 0.1 0.15 0.2 0.25 


(N-Z)/A 


FIG. 4: (Colour online) Same as Figj3l but for K J 4 of spherical A = 100,132 isobar with respect 
to (N — Z)/A. 


C. Sn isotope 


We obtain the expansion coefficients of K& in previous sections. To confirm the coef¬ 
ficients and the expansion Eq.(J7J) itself, we compare the Eqmr in Eq.(| 6 ]) with the coeffi¬ 
cients (K sut{ , K r , Kqoui), in the experiment and Egmr by the linear response calculation 
with Eq. flTOlh for Sn isotopes (H=100 - 132). Figure [5] shows the E GMR for Sn isotopes. 
Solid lines and filled symbols mean the results directly calculated with Eq. flTOlh dashed lines 
correspond to the Eqmr in Eq.(( 6 ]) with the coefficients of case (i) for each interaction, and 
open square symbols are experimental data at RCNPjlOj. Furthermore, to show the pairing 
effects we add the E G mr with TDHF only for SkM*, which are symbolized by open circles 
and dotted line. 

The dashed lines well reproduce whole E'gmr of Sn isotopes within a smaller than 0.3 
MeV. It means the expansion Eq. flT]) is an effective procedure and the coefficients are suitable. 
The comparison our results and he experimental data may recommend SkM* and SGII 
parameters as a candidate of the “answer 3 '. 

The pairing effects to the trend on Sn isotope ISGMR can be discussed in the comparison 
between our results obtained by TDHF (open) and by Cb-TDHFB (filled) with SkM*. The 
small difference between them appears in whole isotopes, which can be expected due to the 
small deformation of HF ground states. While the effect to soften EOS slightly in a surface- 
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type pairing functional was reported 


32 


35 


, our results indicate the opposite effects, 


however whose mechanism is different from the previous studies. The HF ground states in 
Sn isotopes have some deformed aspects as Table [TQ As the explanation in Sec. IIII A1 the 
centroid energy of ISGMR is estimated at lower than that of spherical nucleus due to the 
coupling with other modes, while using the summation analysis such as Eq. lfTUl) . 



FIG. 5: (Colour online) Same as Figj2j but for Sn isotopes from A =100 to 132 with respect to 


A. Filled and open circle means the result calculated by Cb-TDHFB and by TDHF with SkM*, 


respectively. Square symbols shows experimental data at RCNP 


Q. 


IV. DISCUSSION 


We obtain the expansion coefficients (K su r f, K r , Rcoui) through the Sec Jill Al and IIIIB1 
and confirm them in Sn isotopes. SkM* and SGII parameters might be likely candidates in 
Sec J 111 (1 however it is not clear that they can approach to more heavier system, at least in 
present results. The many GMR distributions for whole mass region have been measured 


in the past, which are summarised in Ref. 


34j . We extend the mass region to apply Eq.( 


and attempt to narrow down the parameters. 

Figure El shows exp erimental Fcmrs (square) which are measured at Texas A&M, Gren- 
noble, Groningen |34| and RCNP|10l. Ill|. and the centroid energy (solid line) by Eq.((6]) for 
each Skyrrne parameter but with (r 2 ) expressed as \J(r 2 )= 0.895R 1 / 3 + 0.321 fin which is 
obtained by fitting systematic mean radii calculated in SkM*. The experimental data cover 
a wide mass region A= 24—238 with Z= 12—92. In A <50 mass region, the mass dependence 
of experiment values is not clear apparently. Over <1=100, they make a visible trend, though 
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FIG. 6: (Colour online) The peak positions of GMR in experiment (square) and the lines obtained 
from Eq.([6|) with the Ka for each Skyrrne parameter, where the the mean radii (r 2 ) are calculated 
with SkM* parameter, are shown. 


some values around A=90 deviate from the trend. 

Whole mass region, the solid lines keep their magnitude order which is same one in FigJHJ. 
thus namely SkI3, SkT3, SLy4, SGII and SkM* in decreasing order. For SkI3, the results 
overestimate experiments. For SGII and SkM* they underestimate over A=160, while they 
can well reproduce those around A=110. If we emphasize the agreement with a heavy 
system, SGII and SkM* can not reach to the conclusion as the best choice. SkM should be 
excluded from the candidate, due to the missing over A=100. From the agreements with 
Eqmr over A=140, SLy4 and SkT3 are good candidates. Therefore our results indicate that 
the range of K r is -305±10 MeV, and the incompressibility parameters for infinite system 
corresponding to the K r are A") x> =-340±35, K 00 = 225±11 and iF sym =-138±18 MeV, which 
are expected from SGII, SkM*, SLy4 and SkT3 parameters. 

In this paper, our analysis has constructed within the assumption Eq. (l6]l and (J7J). As 
Figj6] shown that, the expansion Eq.([7|) well work to connect between the infinite nuclear 
properties and finite nuclear excitation mode in Eq.(j6]). However, the mass dependence 
of Agmr in experiment is apparently different from that scaled by A -1 / 3 , which has been 


mentioned also in Ref. 


36 


37]. The small panel in Figj6] shows the fitting experimental 
and by A _1 / 3 + A _1 / 6 (solid). The coefficient of A -1 / 3 


data for A=90—238, by A -1 / 3 (dashed 
is consistent with the ordinary value 
region. The solid line reproduces experimental data in a whole mass region. The A” 1 / 3 and 


30], but. it underestimates experiment in heavy mass 
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A -1 / 6 scaling are well known on the GDR as Steinwedel-Jensen (SJ) and Goldhaber-Teller 
(GT) model |38|, respectively in which the dipole moments are caused by polarization and 
by proton-neutron displacement. If we adapt the models to GMR, the SJ and GT-model 
will respectively correspond to the density wave oscillation keeping nuclear size, and to the 
vibration between expansion and contraction of the size. 


V. CONCLUSION 

We have performed a systematic investigation of the ISGMR to extract the isospin de¬ 
pendent compression EOS parameter from finite nuclear system, using the linear response 
calculation with Cb-TDHFB represented in 3D coordinate space. The expansion coeffi¬ 
cients (Rgurf, K t , K Cou \ ) of finite incompressibility Ka are determined from N=Z and 
isobar A=100, 132 nuclide for each Skyrme interaction (SGII, SkM*, SLy4, SkT3, SkI3). 
Furthermore in Sn isotopes (A=100—132), it is confirmed that the coefficients are suitable 
and the Ka expansion is available. 

The absolute values of K sur { are lager than those of AGo, thus the analysis assuming 
K sm i——K 00 is not effective, which are often used in previous study for non-relativistic 
interaction. The Acoui estimation Eq.([TT|) works well. The isospin term K r in the hnite 
system is smaller than the K%°. The magi city in a present work does not affect the conclusion. 
And the pairing effect in Sn isotopes are not so large, although the deformation due to 
the lack of pairing causes the coupling quadrupole and monopole modes, and disturbs the 
position of GMR. 

We narrow down the values from the comparison experimental data with our results in a 
whole mass region, which indicates the range of K r is -305±10 MeV. The incompressibility 
parameters corresponding to the range have Af X) =225±ll and A' sym =-138±18 MeV, which is 
consistent with the previous study of slop parameter L As an indication of Eq. (fT2]h the 
L has an important role to decide Rs ym . Thus the studies related to L should be also pressed 
forward in parallel. The ordinary mass scaling of GMR has underestimation in heavy mass 
system, which may indicate the needs to reconsider the expression of incompressibility. 

To determine the isospin term AT r , the systematic ISGMR data in isobar nuclide are very 
effective. They should be heavy system over A=100 at least, because in light system the 
nuclear speciality is showing up. When the heavy nucleus is used to analyse GMR, we should 
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note the deformation of open shell nuclei. If the monopole and other modes in deformed 
nucleus can be separated, the analysis data increases and the relation EOS parameter and 
finite system will be more robust. 
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TABLE I: Root mean square radius R = \J(r 2 ) [fm] and absolute value of quadrupole deformation 
parameter f5 in ground state for N = Z and A = 100,132 nuclide. 



SGII 

SkM* 

SLy4 

SkT3 

SkI3 


R |/3 

R |/3| 

R |/3 

R \(3\ 

R \(3\ 

40 Ca 

3.35 0.00 

3.40 0.00 

3.39 0.00 

3.40 0.00 

3.36 0.00 

44 Ti 

3.44 0.00 

3.48 0.00 

3.48 0.00 

3.47 0.00 

3.46 0.00 

48 Cr 

3.56 0.26 

3.58 0.21 

3.60 0.26 

3.53 0.03 

3.58 0.29 

52 Fe 

3.58 0.00 

3.61 0.00 

3.62 0.00 

3.60 0.00 

3.60 0.00 

56 Ni 

3.64 0.00 

3.66 0.00 

3.67 0.00 

3.65 0.00 

3.65 0.00 

60 Zn 

3.77 0.20 

3.75 0.00 

3.77 0.02 

3.74 0.00 

3.76 0.02 

64 Ge 

3.83 0.00 

3.84 0.00 

3.86 0.00 

3.82 0.00 

3.86 0.00 

68 Se 

3.95 0.22 

3.94 0.09 

3.97 0.17 

3.91 0.01 

3.97 0.23 

72 Kr 

4.02 0.19 

4.07 0.26 

4.05 0.17 

4.03 0.23 

4.08 0.31 

76 Sr 

4.06 0.01 

4.12 0.00 

4.08 0.01 

4.15 0.39 

4.08 0.13 

80 Zr 

4.12 0.00 

4.15 0.00 

4.15 0.00 

4.12 0.00 

4.12 0.00 

84 Mo 

4.17 0.00 

4.20 0.00 

4.20 0.00 

4.17 0.00 

4.18 0.00 

00 

00 

S0 

c 

4.22 0.00 

4.25 0.00 

4.25 0.00 

4.23 0.00 

4.23 0.00 

92 Pd 

4.27 0.00 

4.29 0.00 

4.30 0.00 

4.27 0.00 

4.28 0.00 

96 Cd 

4.32 0.00 

4.34 0.00 

4.34 0.00 

4.32 0.00 

4.32 0.00 

100 Sn 

4.36 0.00 

4.38 0.00 

4.39 0.00 

4.36 0.00 

4.36 0.00 

100 Kr 4.50 0.23 

4.54 0.23 

4.54 0.22 

4.50 0.21 

4.57 0.24 

100 Sr 4.54 0.39 

4.56 0.38 

4.57 0.39 

4.52 0.36 

4.59 0.40 

100 Zr 4.52 0.38 

4.52 0.36 

4.54 0.36 

4.48 0.33 

4.59 0.43 

100 Mo 4.40 0.00 

4.42 0.00 

4.43 0.00 

4.38 0.00 

4.46 0.20 

100 Ru 4.39 0.00 

4.40 0.00 

4.42 0.00 

4.37 0.00 

4.43 0.19 

100 Pd 4.38 0.00 

4.39 0.00 

4.40 0.00 

4.36 0.00 

4.39 0.01 

100 Cd 4.37 0.00 

4.38 0.00 

4.39 0.00 

4.35 0.00 

4.38 0.00 

132 Sn 4.78 0.00 

4.80 0.00 

4.80 0.00 

4.79 0.00 

4.82 0.00 

132 Te 

4.79 0.00 

4.80 0.00 

4.81 0.00 

4.78 0.00 

4.82 0.00 

132 Xe 

4.79 0.00 

4.81 0.00 

4.81 0.00 

4.78 0.00 

4.82 0.00 

132 Ba 4.81 0.15 

4.83 0.161^.84 0.15 

4.79 0.14 

4.83 0.16 


132 Ce 4.83 0.21 4.85 0.23 4.86 0.21 4.82 0.23 4.84 0.21 













TABLE II: Same as Table Q] but for Sn isotopes with from N =52 to 80. Although the HF+BCS 
ground states take only a spherical shape in present work, the HF results in SkM* have some 
deformation. 



SGII 

R 

R 

SkM* 

7?hf 

|/3hf 

SLy4 SkT3 SkI3 

R R R 

102 Sn 

4.39 

4.405 

4.401 

0.04 

4.42 

4.38 

4.40 

104 Sn 

4.43 

4.434 

4.426 

0.05 

4.45 

4.40 

4.43 

106 Sn 

4.46 

4.462 

4.448 

0.00 

4.48 

4.43 

4.47 

108 Sn 

4.48 

4.491 

4.488 

0.07 

4.51 

4.46 

4.50 

110 Sn 

4.51 

4.521 

4.527 

0.10 

4.54 

4.49 

4.53 

112 Sn 

4.54 

4.551 

4.558 

0.10 

4.57 

4.52 

4.56 

114 Sn 

4.57 

4.581 

4.596 

0.11 

4.60 

4.54 

4.59 

116 Sn 

4.59 

4.609 

4.632 

0.15 

4.62 

4.57 

4.61 

118 Sn 

4.62 

4.636 

4.699 

0.28 

4.65 

4.60 

4.64 

120 Sn 

4.64 

4.662 

4.677 

0.17 

4.67 

4.63 

4.67 

122 Sn 

4.67 

4.688 

4.703 

0.14 

4.70 

4.66 

4.70 

124 Sn 

4.69 

4.712 

4.724 

0.10 

4.72 

4.68 

4.72 

126 Sn 

4.71 

4.736 

4.743 

0.00 

4.74 

4.71 

4.75 

128 Sn 

4.74 

4.759 

4.766 

0.06 

4.76 

4.74 

4.77 

130 Sn 

4.76 

4.781 

4.784 

0.04 

4.79 

4.76 

4.80 
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TABLE III: Parameters related to EOS (po[fna~ 3 ], Qo, L, A sym , K%°, Kcoul [MeV]) and finite 
incompressibility (AT sur f, K T , ILcoul [MeV]) for each Skyrme interaction. The K^° and Acoul are ob¬ 
tained by Eq. (fT2]h(fTT]h Nucleon mass mc 2 =938.9187 MeV, hc=197.327 MeV frn and a _1 =137.036 


are used 


31]- 


EDF 


(i) without DM (ii) with DM 


Int. po K 00 Q o L K sym K T K-Cou i 
SGII .1583 214.6 -380.9 37.63 -145.9 -304.9 -4.69 
SkM* .1603 216.6 -386.1 45.78 -155.9 -349.0 -4.70 
SLy4 .1595 229.9 -363.1 45.96 -119.7 -322.9 -4.85 
SkT3 .1610 235.7 -382.7 55.31 -132.1 -374.1 -4.83 
SkI3 .1598 258.2 -303.8 100.5 73.03 -411.8 -5.13 


A-surf h A Coul 

-273.7 -295.9 -4.52 
-285.3 -282.0 -4.76 
-318.4 -314.3 -4.73 
-326.1 -305.1 -4.54 
-372.6 -352.0 -4.75 


Asurf ^Ooul 

-260.4 -283.5 -5.10 
-276.7 -285.0 -5.09 
-308.3 -305.4 -5.13 
-325.0 -305.2 -4.51 
-364.8 -332.6 -4.99 
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